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ABSTRACT 

By combining data from the NASA Wide-field Infrared Survey Explorer ( WISE) mission with optical spectroscopy 
from the W. M. Keck telescope, we discover a mid-IR color criterion that yields a 78% success rate in identifying 
rare, typically radio-quiet, 1.6 < z < 4.6 dusty Lya emitters (LAEs). Of these, at least 37% have emission 
extended on scales of 30-100 kpc and are considered Lya “blobs” (LABs), The objects have a surface density of 
only ^0, 1 deg” 2 , making them rare enough that they have been largely missed in deep, small area surveys. We 
measured spectroscopic redshifts for 92 of these galaxies, and find that the LAEs (LABs) have a median redshift 
of 2.3 (2.5). The WISE photometry coupled with data from Herschel ( Herschel is an ESA space observatory with 
science instruments provided by European-led Principal Investigator consortia and with important participation 
from NASA) reveals that these galaxies are in the Hyper Luminous IR galaxy regime (Ljr > 10 13 -10 14 L©) and 
have warm colors. They are typically more luminous and warmer than other dusty, z ~ 2 populations such as 
submillimeter-selected galaxies and dust-obscured galaxies. These traits are commonly associated with the dust 
being illuminated by intense active galactic nucleus activity. We hypothesize that the combination of spatially 
extended Lya, large amounts of warm IR-luminous dust, and rarity (implying a short-lived phase) can be explained 
if the galaxies are undergoing brief, intense “feedback” transforming them from an extreme dusty starburst/QSO 
into a mature galaxy. 
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1. INTRODUCTION 

High-redshift Lya emission is widely used to study star 

formation in galaxies (e.g., Cowie & Hu 1998; Steidel et al. 

2000; Gawiser et al, 2007; Gronwall et al. 2007; Nilsson et al. 

2007; Finkelstein et al. 2009; Ono et al. 2010). Systems that 

exhibit this line are typically referred to as Lya emitters (LAEs); 

a small subset show extended emission on scales >30 kpc 

(occasionally upward of 100 kpc) and are considered Lya 

“blobs” (LABs). Among the largest coherent galactic structures 

known in the universe, LABs are extremely energetic, with 

Lya luminosities of ~10 42 — 10 44 erg s’ 1 , and have been shown 

to trace over-densities of galaxies at high-redshift (e.g.. Keel 
et al. 1999; Francis et al. 2001; Steidel et al. 2000; Matsuda 
et al. 2004; Prescott et al. 2008, 2012a; Yang et al. 2009; Erb et al. 
201 1). Bright, optically selected LABs are 100-1000 times less 
abundant than LAEs (e.g., Steidel et al. 2000; Matsuda et al. 
2004; Saito et al. 2006) with only a few dozen known (Steidel 
et al. 2000; Matsuda et al. 2004, 2011; Dey et al. 2005; Geach 
et al. 2005; Nilsson et al. 2006; Lai et al. 2007; Smith & Jarvis 

2007; Prescott et al. 2009, 2012b; Yang et al, 2009, 2010; Ouchi 
et al. 2009; Webb et al. 2009). 

Resonant scattering by atomic hydrogen means Lya photons 
are easily extinguished by dust. This suggests that dusty infrared 


galaxies and LAE/LAB populations, despite both peaking in 
activity around z ^ 2, should largely not overlap. Indeed Spitzer 
observations of ~40 optically selected z = 2-3.5 LABs found 
that only <15% are also luminous in the infrared (Webb et al. 
2009). Furthermore, no known optically discovered LAB is 
WISE detected at 12 or 22 /zm (see Section 3.1 of this work) 
including the first LAB found through its dust emission at 24 pm 
(Dey et al. 2005). Although, 40%-50% of the infrared luminous 
submillimeter galaxies (SMGs) show Lya emission, they rarely 
exhibit spatially extended Lya emission (Chapman et al. 2005). 
This suggests however, that dust and Lya emission are not in 
fact mutually exclusive, and points to a clumpy dust distribution 
that allows sight lines for Lya photons to escape (Neufeld 1991). 
Since SMGs and related high-redshift galaxies are thought to 
be just one stage in the evolution of massive elliptical galaxies, 
the study of spatially extended Lya in dusty galaxies can thus 
provide a unique insight into this process. 

This discovery paper is the first in a series, and presents a new 
and efficient mid-IR color technique to select a never before seen 
population of z ^ 2-4 dusty LAEs, of which at least one-third 
are also found to be LABs. We use the redshift distribution, 
general spectroscopic properties, and preliminary Herschel 
observations to determine how these objects relate toother dusty 
galaxies at these redshifts and speculate how these enormous 
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Figure L Two- dimensional LR1S spectra for four of the WLABs. The spatial exieni of ihe I emission is noted in the corner of each image. These examples are 
representative of the asymmetric, spatially extended line profile, and clumpy nature of the Lyn emission found in these galaxies* The arrows are meant to guide the 
eye and highlight the variation in velocity profiles. 

(A color version of this figure is available in the online journal) 


LABs are being powered. We assume a ACDM cosmology with 
Q m = 0,27* 0 A = 0.73, and Ho = 71 kin s _l Mpc -1 and use 
Vega magnitudes, 

2. SAMPLE DEFINITION AND OBSERVATIONS 

One of the primary science goals of the NASA Wide -fie Id 
Infrared Survey Explorer (WISE) mission (Wright et al. 2010). 
which covers the whole sky at wavelengths of 3,4, 4,6, 12, and 
22 /im (W1-W4), is to find the most luminous galaxies in the 
universe. We performed initial searches in WISE color space 
which revealed a diversity of galaxy types over a wide range 
of redshifts (Griffith et al, 20 M: Eisenhardt et al. 2012: Wu 
et al, 2012: Lonsdale et ah 2013; C, R. Bridge et al. 2013* 
in preparation). In this discovery' paper we report the results 
of a deep spectroscopic follow-up campaign to target optically 
faint galaxies with the reddest WISE (mid-IR) colors. A striking 
finding of this work is that not only did the majority of these 
high-redshift sources have intense Lyot emission at least one- 
third of these U7SE-selected LAEs are LABs, having spatially 
extended Lya reaching 30-100 kpc (see Figure 1). 


We refer to this new population of hot, dusty, W/S£-detected 
LAEs/LABs as — WLAEs and WLABs because they are se- 
lected by WISE colors, and cover a previously unexplored wave- 
length and flux density parameter space over the whole sky. 

2, /, Selection Method 

Realizing that WLAEs and WLABs lie in a very specific 
region of a WISE color-magnitude diagram, we refined the 
selection criteria used in the initial spectroscopic follow-up to 
focus on high-redshift LAEs/LABs. 

The resulting ”WLAB” selection criteria, shown in Figure 2, 
requires W2-W3 (4.5-12 ftm) ^ 4,8 and S/N > 5 in W3 and 
W4. We also excluded sources within 30 of the Galactic center 
and 10 of the Galactic plane to avoid asymptotic giant branch 
stars and saturation artifacts. In order to remove low-redshitt 
(z < 0.5) star- forming galaxies with similar H7S£ colors, 
we required non-detections in Sloan Digital Sky Survey and 
Digitized Sky Survey imaging ( r f > 22). Finally, each object's 
co- added W3 and W4 images were visually inspected to ensure 
the source is not spurious. These criteria result in an average 
source density of ^0. 1 deg - ^. 
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Figure 2, W7S£ W3(12 /cm) vs. WISE W2-W3 {4.5-12 /cm) color for all 
sources in the sample with spectroscopic redshifis between 1.6 ^ z K 4,6, 
WLAEs (black circles), 30-40 kpc WLAUs (green squares), and 40- 1 (X) kpe 
WLABs (orange squares) are plotted. For comparison. z = 2-3 WISE detected 
SMCis/DOGs ;ire highlighted: Magnelli el ah (20 L 2* open squares), Bussmann 
et al. (2012. cyan downward triangle). Chapman et ah (2010. black asterisks). 
Smith et ah (2009, gray triangle). Chapman et id. (2005. red circle), and 
Chapman et ah (2005, blue star). LABdOS is not detected by WISE, hence 
the data shown is from Spirzer (Dey el ah 2005). The WLAEs (WLABs) exhibit 
unique, redder mid-IR colors (implying hotter dust temperatures) than other 
dusty z ^ 2 populations. The WLAE/WLAB selection technique has a 78% 
(>29%) success rate in identifying ; > 1.6 dusty IAEs (LABs), 

{A color version of this figure is available in the online journal.) 


2.2. Optical Spectroscopy 

To determine redshifts and discern the nature of the galaxies, 
we used Keek I LR1S (Oke et ah. 1995) over the course of six 
runs between 2010 July and 2012 January'. The sample contains 
101 objects that fulfill die WLAE/WLAB selection, of which 
dose to half had previously been followed-up using other red 
selections (i.e., \V1 W2-drop criteria; Eisenhardl el ah 2012) and 
therefore did not need to be re-observed. The full spectroscopic 
campaign will be presented in an upcoming paper. The majority 
of data were taken using the 600 i mm -1 grism in the blue arm 
(Abira — 4000 A; spectral resolving power R = a/Aa ^ 750), 
the 400 t mm” 1 grating on the red arm (k b \ :ui > = 7800 A: 
R — 700), the 5600 A dichroic and a l '.'5 wide long slit. The 
targets were observed with a median seeing of O'' 7 with total 
integration times between 2CW0 minutes. The LR1S data were 
reduced and flux calibrated using standard procedures. 

The optical spectroscopy reveals a wide range of character- 
istics, from pure starburst galaxies with narrow emission lines 
and several ultraviolet interstellar absorption lines comparable 
to Lyman break galaxies (Steidel et ah 2000) to those with strong 
active galactic nucleus (AGN) components, inferred from the 
presence of high ionization lines such as N v, Si iv/O iv T C iv 
and broad emission line profiles of several thousand km s”L 

Robust redshifts based on two or more spectral features 
(primarily N v, Siiv/Oiv, Civ, Mg n, Nett) were determined 
for 92 out of ihe 101 sources we targeted (91% success rate), 
while the remaining 9 showed no continuum or spectral lines. We 
believe these nine sources are most likely highly dust-obscured 
and hence below our detection threshold. The sample has a 
redshift range of z = 1,13-4.59 (see Figure 3). The WISE 
LAEs (LABs) have a roughly Gaussian redshift distribution 
centered around c — 2.5 (23), The excellent blue sensitivity 
of LR1S allows Lya emission to be detected down to 3160 A, 
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Figure 3. Spectroscopic redshift distribution for ihe full sample 
(WLAHs+WLABs; tilled gray), and WLABs (filled orange). For compari- 
son we plot the redshift distribution of other well-studied dusty populations: 
SMGs (black dashed: Chapman et at, 2005), Hmchet SMGs ( blue long dashed: 
Magnelli et al. 2012), and DOGs (red dotted; Bussmann et al. 2012), The VV7SF 
LAEs (LABs) have a roughly Gaussian redshift distribution centered around 
Z = 2.5 (2.3), very' si milar to that of SMGs and DOGs, 

{A color version of ibis figure is available in the online journal.) 


corresponding to a redshift of z ^ 1 ,6, Of the 92 lV7S£-selecied 
galaxies with a robust redshift, 3 were c < L6 star-forming 
galaxies, while 89 were at a c ^ L6, and of those 79 (89%) 
showed Lycr in emission. Given all of the spectra exhibiied 
multiple emission lines (typically Lya, N v, Civ), a secure 
idemification of Ly a was possible. An additionally striking 
finding was lhar 37% (29/79) of the WLAEs showed Lye* 
emission in the two-dimensional spectra extended on spatial 
scales of ^30 kpc, with 1 8% (14) having emission 40-100 kpc 
in extent (i.e,, LABs). 

Given we used random slit orientations the detection rate is 
consistent with all the LAEs in our sample being are LABs. We 
also note that an additional 14 (of the 79 WLAEs) have extended 
emission of 25-30 kpc, which is more than twice the seeing and 
4-5 times the typical sizes of dusty z ^ 2 galaxies (i.e., SMGs 
and dust-obscured galaxies (DOGs); Chapman et al, 2005; Biggs 
& Ivison 2008; Bussmann et al. 2009; Casey et al. 2009; 
lono et al. 2009; Bothwell et al. 2010; Younger et al, 2010), 
These therefore constitute smaller LABs or Ly a clouds. 

The definition of a LAB does vary between studies. 
Matsu da et al. (2004) required an isophotal area of at 
least 16 arcsec 2 (or —30 kpc in length at z = 3,1) 
and a surface brightness of ^27,6 mag arcsec" 2 {^2 x 
|Q“ 18 erg s" 1 cm -2 arcsec -2 ), while Matsuda et al, (2011) 
probed down to 1,4 x 10 -IS erg s -1 cm” 2 arcsec -2 . In contrast to 
these narrow-band imaging studies we used the one-dimensional 
spectra to measure the spatial extent of the Lyo emission in one 
direction. Multiple slit orientations were used to probe the spa- 
tial line morphology for four of the LABs. A slit position angle 
rotated by 45 and 90° of the original angle often showed little or 
no extended emission, implying an asymmetric and filamentary' 
line morphology. 

The depth of the spectroscopy reached Lya surface bright- 
nesses of 1-100 x 10“ 18 erg s” 1 cm” 2 aresec” 2 . This is well 
within the range used in optical narrow-band surveys to iden- 
tify LABs (e,g., Steidel et al. 2000; Matsuda et al. 2004, 2011; 
Nilsson et al. 2006; Smith Sc Jarvis 2007; Ouchi et al. 2009; 
Yang et al, 2009* 2010). Each WLAE spectra was inspected and 
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Table 1 

Far- IK Properties of the Eight Objects with Herschel Observations 


z 

5(70) 

5(160) 

5(250) 

5(350) 

5(500) 

T 

(K) 

log(t/L 0 ) 

2.022 

27 

80 

63 

52 

43 

52 ± 8 

13.1 ±0.1 

2.140 

84 

217 

172 

110 

53 

70 ± 6 

13,8 ±0.1 

2.298 

40 

75 

109 

82 

44 

46 ± 6 

13,3 ±0,1 

2.450 

18 

60 

62 

41 

<30 

60:1:7 

13.2 ±0.1 

2.452 

29 

53 

47 

33 

<30 

84 ±8 

13.4 ±0.1 

2.542 

24 

40 

39 

34 

<30 

84 ± 12 

13.3 ±0.1 

3.048 

8 

37 

59 

56 

43 

48 ±3 

13.3 ±0.1 

3.228 

16 

56 

53 

45 

<30 

71 ± 10 

13.5 ± 0.1 

Notes. 

Flux errors (la) 

lor the 

Herschel 

bands are 4. 10. 

6. 6, 8 mJy 


corresponding u> 7(1 160, 250, 350. and 500 Dust temperature and 
luminosity were lit using a gray body with 5{i') ct Biv)v l where £t(M) is 
the Planck function and the statistical la errors are quoted, 

those with spatially extended Lyot emission >30 kpc were clas- 
sified as WISE LABs or WLABs. Given that the Lya surface 
brightness reached by the spectra conservatively satisfy, and in 
many cases are brighter than the definitions of traditional LABs, 
it is highly plausible that the mie sizes of the WLABs are in fact 
larger under these standard definitions. The sizes of the individ- 
ual WLABs and their line fluxes will be presented in depth in 
the catalog paper (C. R. Bridge et al. 2013, in preparation). Ap- 
plying the LAE (LABs) fractions found by our spectroscopic 
survey to the total number of sources that satisfy the color 
criteria we estimated a source density for WLABs (WLABs) 
of M). 1 (0,03) deg -2 , or approximately 4500 ( 1600) over the 
full sky. 

23, Far- infrared Photometry 

The red WISE colors strongly argue that significant amounts 
of warm dust are present in these galaxies. Since the spectral 
energy distributions (SEDs) of dusty galaxies typically peak in 
the far- infrared, observations between 0. 1-1 nun are required in 
order to accurately determine the total IR luminosity and dust 
temperature of the sources. To that end, we have undertaken 
a campaign with Herschel to study the WLABs in five bands 
covering 70-500 /xm using the FACS (Poglitsch et al. 2010) and 
SPIRE (Griffin et al, 2010) imagers. 

Eight objects with WISE colors drawn randomly from our 
selection have been observed with Herschel to date, and all 
were detected in at least four bands (see Table l). To estimate 
the dust luminosity and temperature, these data were fit using a 
modified blackbody with S(v) a B( i\ T)3 Here, B(v , T) is 
the Planck function. 

While more complex models can be used, this parameteriza- 
tion has been adopted by many other studies, enabling a straight- 
forward comparison with, for example, samples of SMGs and 
DOGs presented in Magnelli et al, (2012) and Melbourne et al, 
(2012), respectively. Figure 4 presents the results of these fits 
and indicates that WLABs typically contain warmer dusi than 
ihese other samples. This suggests that they are being heated, at 
least in part, by an AGN, since starburst galaxies such as SMGs 
are typically colder (also see Wu et al. 2012). The WLABs 
are systematically more luminous than SMGs, with bolomet- 
ric luminosities in excess of Lrk > !0 13 Lq* making them 
Hyper- Luminous Infrared Galaxies (HyLIRGs). 

It is tempting to associate these high luminosities with 
gravitational lensing; however, none of the optical spectra 
demonstrate contamination from a foreground source, and the 



Figure 4, Estimated dust properties based on Herschel observations. Objects 
between 2,0 ^ z ^ 3.5 extracted from catalogs of DOGs and unlensed SMGs 
are plotted against the eight observed WLABs which cover a similar redshift 
range. In all cases, fits are to a single temperature modified blackbody, with a 
fixed emissivhy of ft = 1.5. Error bars in all cases are based on fit results only 
and do not account for oilier systematic effects. 

(A color version of this figure is available in the online journal.) 

spatially extended Lya emission for all the sources is at the 
same redshift as the oilier spectral features. Furthermore, visual 
inspection of near- 1 R, adaptive optic (AO) imaging with Keck's 
NIRC2 camera (Eisenhardi et al. 2012; C. R. Bridge et al. 
2013, in preparation), and Hubble Space Telescope WFC3/JR 
imaging of over a dozen WLABs (C. R. Bridge ei ai. 2013, in 
preparation) reveal no evidence of either a potential lens (nearby 
galaxies or clusters) or any morphological signatures, such as 
arcs or multiple components that the WISE object itself is lensed. 
However a full lensing analysis is required to rule out any small 
amplifications. 

3. DISCUSSION 

3.1. WISE LABs and Previously Known LABs 

One of the main features that sets this new population of 
WLAEs apart from other high-redshift galaxies, aside from 
their extremely red mid-lR color, is the unprecedented fraction 
that exhibit spatially extended Lya emission, despite having 
large amounts of dust. Traditionally, LABs have been discov- 
ered sere n dip itou sly through narrow-band optical surveys that 
concentrate on a particular redshift (e.g., Steidel et at. 2()00; 
Maisuda et al. 2004, 2011; Nilsson et al. 2006; Smith & Jarvis 
2007; Yang et al. 2009, 2010; Ouchi et al. 2009), and more re- 
cently using optical broad-band surveys (Prescott et al. 2012b, 
2013). These types of surveys more closely resemble “blind” 
searches, in that there is no a priori knowledge of where or how 
many LABs reside in the field. In contrast, the technique pre- 
sented here {Figure 2) uses broad-band photometric data alone 
to preselect with high efficiency a population of LABs over a 
wider range of redshifts, with the caveat that they are seemingly 
much dustier than the ones found in optical surveys. 

Although WLABs. at the correct redshift, would be de- 
tected in most optical narrow-band LAB searches, this type 
of serendipitous methodology can only efficiently patrol a small 
area, from a few square aremins up to a couple of square de- 
grees of the sky (Steidel et al. 2000; Matsuda et al. 2004, 2011). 
Given a WLAB surface density of 0.03 deg -2 , the WLABs 
are thousands of times rarer than their optical LAB cousins. 
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meaning that we should not have expected to discover any 
WLABs in existing optical LABs surveys. Again, this is con- 
sistent with our finding, that no optical LAB has the extreme 
mid-lR colors exhibited by the WLABs. We proceed now with 
the details of that investigation. 

We cross matched the WISE catalog with all previously 
known LABs in the literature (Steidel et al. 2000; Matsuda et al. 
2004, 2011; Dey et al. 2005; Geach et al. 2005; Nilsson et ah 
2006; Lai et al. 2007; Smith & Jarvis 2007; Prescott et al. 2009, 
2012b; Yang et al. 2009, 2010; Ouchi et al. 2009; Smith et al. 
2009; Webb et al. 2009) and found that none are detected at 
12 or 22 /tm above a S/N of 5. It is important to note that 
this includes all published LABs known to have mid- or far- 
IR detections from follow-up observations with Spitzer and 
Herschel (Geach et al. 2005; Dey et al. 2005; Webb et al. 2009). 
No previously known LABs, including those which exhibit IR 
mission, are WISE sources, and would not be selected by the 
WL AB color technique. It is clear that WLABs occupy a distinct 
region of mid-IR color space, and are likely in a different stage 
of evolution, environment, or energy generation process (i.e., 
AGN, extreme starburst) then previously discovered LABs. 

Despite the WLABs being a factor of 10-1000 times more 
luminous in the mid-far-IR then any previously discovered 
LAB, they have Lya luminosities (10 4 -10 44 erg s' 1 ), rest- 
frame equivalent widths (a few tens to 300 A), and LABs sizes 
similar to those discovered optically. In other words, in the 
optical WLABs are similar to LABs but they have a larger 
fraction of mid-IR power to Lya emission. We hypothesize that 
the success of the WLAB technique is due to it being able to 
select AGN at or near its peak of energy output. Such an object 
could produce strong outflows clearing paths through the dust 
for the Lya photons to escape. Indeed our optical spectra do 
exhibit lines consistent with a strong outflow (see Section 3.3 
for a discussion). While supemovae (SNe) are also thought to 
provide a mechanism for clearing out dust, a powerful starburst 
alone cannot account for the strong mid-IR luminosities, red 
colors or hot (50-85 K) dust temperatures (Section 2.3 and 
Table 1 ) displayed by these systems. Additionally, SMGs, which 
demonstrate similar levels of star formation do not exhibit 
extended Lya emission. Therefore, intense star formation is 
likely not the primary process responsible for the WLABs. 

It is worth reiterating that the phenomena of large Lya halos 
has also been found around some high-redshift radio galaxies 
(HzRGs) and QSOs (e.g., McCarthy et al. 1987; Heckman 
et al. 1991; Reuland et al. 2003; Smith et al. 2009). In these 
systems the extended Lya emission is thought to be powered by 
interactions between the radio jets and the ambient intergalactic 
medium (IGM). The majority (two-thirds) of our sources are 
undetected in the FIRST (< 1 mJy) or NVSS radio surveys, and 
while one-third do show faint radio fluxes they are well off the 
Lya-radio correlation (in both Lya and radio flux) found for 
HzRGs (Smith et al. 2009), differentiating them from this well- 
studied population. However it may be the case that the WLABs 
are in a radio quiet stage and may evolve into similar systems. 

3.2. WLABs and Other Dusty High-redshift Populations 

Though the IR luminosity and redshift distribution of 
the WLAEs/WLABs are similar to other populations of 
IR-luminous galaxies (Figure 3), they must differ in some way 
given the large fraction with extended Lya emission. 

Blind submillimeter surveys have been very successful 
at finding high-redshift, dust-obscured galaxies with far-IR 


luminosities in excess of 10 12 L Q (i.e., SMGs; e.g., Smail 
et al. 1997; Barger et al. 1999; Blain et al. 1999; Borys et al. 
2005; Hainline et al. 2009). However, these long wavelength 
observations (>500 /zm) select systems with cooler dust and 
less extreme mid-IR colors than seen in our WISE sample 
(Figures 2 and 4). Generally, SMGs are powered by star for- 
mation rather than Compton-thin AGN (Alexander et al. 2005). 
Based on the limited Herschel observations conducted so far, 
most of the WLABs would also be classified as SMGs. However 
only \%-2% of known SMGs fulfill the WLAE/WLAB color 
criteria (this work; Figure 2), and there are only a couple of rare 
SMGs with extended Lya emission (Chapman et al. 2005). 

Dey et al. (2008) have successfully used a mid-IR to optical 
flux ratio (F^/zm/J 7 /? >1000 and F 24 /Am > 0.3 mJy) to 
select a population of z ~ 2 dust-obscured galaxies, which 
they refer to as DOGs (Houck et al. 2005; Desai et al. 2009; 
Bussmann et al. 2012). By probing blueward of the SED peak, 
this technique recovers a mix of AGN and starburst dominated 
galaxies the majority of which reside at z ^ 2. Given that 
DOGs share a common epoch, dusty composite nature and mid- 
IR selection it makes sense to briefly compare them to the 
WLAEs/WLABs. The DOG sample from Bussmann et al. 
(20 1 2) was crossed checked with the WISE catalog and revealed 
that only -8% of the DOGs (N = 8) fulfill the WLAE/WLAB 
selection presented here. Although there is some overlap in the 
W2-W3 color space. Figure 2 highlights that Spitzer DOGs are 
typically 1.5-2 mag fainter at 12 jim (W3) than the WLAEs/ 
WLABs. Deep optical spectroscopy of ^60% of these Spitzer 
DOGs has been conducted by Bussmann et al. (2012) and one 
source has been reported to exhibited extended Lya (Dey et al. 
2005), however it is not detected by WISE and is more akin 
in mid-IR flux to a traditional SMG. The fainter W3 flux of 
Spitzer DOGs could reflect a lower intrinsic luminosity (also 
see Figure 4) or a different nature. Using the observed r'-W4 
color of the WLAEs/WLABs, we also consider the fraction 
that would satisfy the DOG criterion. All our WLAEs/WLABs 
would be classed as DOGs, although they are too rare to have 
been discovered in Spitzer MIPS 24 surveys. Finally, it is 
important to note that the DOG criterion focuses purely on 
the optical-to-mid-IR flux ratio, and does not specifically select 
sources with a steeply rising mid-IR SED (Figure 2), which 
is what we have shown are generally the objects that exhibit 
extreme Lya properties. 

It is also important to discuss the similarities and differences 
between the WLAE/WLAB selection and other H7S£-based 
color selections that have been undertaken. One such example 
is the “WlW2-dropout$,” outlined in Eisenhardt et al. (2012), 
which has been successful in identifying ULIRGs over a wider 
redshift range (z 0.05-4.6). There is a substantial overlap 
in the z > 2 samples; however, overall we find that the 
refined WLAE/WLAB criterion is twice as effective at selecting 
WLAE/WLABs at high-z. The primary differences between 
these two selections are that the WLAE/WLAB criteria impose 
an optical magnitude limit to remove contamination from low 
redshift galaxies, places no flux cut on W1 (many WLABs are 
detected in W1 and/or W2), probe a factor of 1.2 (2) deeper 
in W3 (W4) flux and are 0.5 mag bluer in W2-W3 than the 
W1 W2-dropout selection. 

Another WISE extragalactic color selection presented in 
Lonsdale et al. (2013), combines specific WISE colors, with 
the requirement of an intermediate or loud radio detection in 
the regions of sky covered by the NVSS 20 cm radio survey 
(Condon et al. 1998). Although these sources share similarly 
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red mid-IR colors to that of the WLAE/WLAB population the 
redshift range of these red radio sources is systematically lower 
(z — 0.4-2.5; Lonsdale et al. 2013) and they show no evidence 
of extended Lyot. As suggested in Lonsdale et al. (2013), 
these red-radio loud objects are good candidates for being in a 
QSO-dominated phase. It is possible that some of the powerful 
AGN thought to reside in both the WLABs and these red- 
radio loud sources are merely at different stages of their 
evolution. We will discuss the potential connection between 
these and the other populations mentioned above in more detail 
in the WLAB catalog paper (C. R. Bridge et al. 2013, in 
preparation). 

3.3. Spectral Properties of WISE LABs 

The cause of the LAB phenomena in optically selected 
samples is still an outstanding question (e.g., Nilsson et al. 
2006; Prescott et al. 2009, 2012a; Yang et al. 2009, 2010; 
Colbert et al. 201 1). A logical explanation is photo-ionization 
by an embedded source, however Matsuda et al. (2004) showed 
that many LABs do not house a sufficiently UV-bright galaxy. 
Optically discovered LABs, like the WISE LABs are radio- 
quiet and hence are not powered by interactions between radio- 
jets and the ambient IGM as in the case of the Lya halos 
surrounding HzRGs (e.g., McCarthy et al. 1987; Heckman 
et al. 1991; Reuland et al. 2003; Smith et al. 2009). Alternative 
mechanisms are ( 1 ) cold accretion, where the LABs are powered 
by gravitational heating during the collapse of a large primordial 
cloud (Haiman et al. 2000; KereS et al. 2005; Dijkstra et al. 
2006; Dekel & Bimboim 2006; Nilsson et al. 2006; Smith 
& Jarvis 2007; Smith et al. 2008), (2) intense starbursts or 
and embedded AGN (e.g.. Keel et al. 1999), (3) shocks from 
SN winds or AGN outflows (e.g., Geach et al. 2005; Colbert 
et al. 2011), and (4) resonant scattering of Lya photons (e.g., 
Steidel et al. 2011; Saito et al. 2006). Given that all the 
WLABs are associated with a luminous IR galaxy and less 
than 15% of optical LABs are (Webb et al. 2009), suggests that 
the dominant powering mechanisms behind these two species 
of LABs are likely different. For clues we investigate the 
optical spectra. 

The high-ionization and broad emission lines seen in >70% 
of optical spectra suggest that the majority, if not all, of the 
WLAEs/WLABs host a dust-obscured AGN, which is also 
consistent with the Herschel and WISE colors. The Lya emission 
is spatially asymmetric and clumpy, and is often offset from 
the systemic redshift of the galaxy by 100-1000’s of km s“ ! , 
implying large-scale outflows. Furthermore, many galaxies in 
the sample demonstrate velocity structure in the Lya line, 
redshifted up to several thousand km s - 1 as a function of distance 
from the central component (Figure 1 (d)). The line profiles are 
diverse, ranging from traditional P-Cygni profiles characteristic 
of Wolf-Rayet and O-star winds to Lya emission peaking on 
the blue side of the Lya profile, suggesting possible inflows 
(Dijkstra & Loeb 2009; Barnes et al. 201 1). 

The filamentary morphology of the extended Lya seen when 
multiple slit orientations were used suggests that the quoted 
efficiency of the color selection in identifying LABs is a 
lower limit. Indeed, using a simple model of the extended 
emission size, geometry, and random slit orientation the 37% 
detection rate is consistent with all the LAEs in our sample 
being LABs. 

The diversity of spectral properties can be explained by either 
varying levels of dust extinction or perhaps the systems are at 
slightly different stages of evolution. 


3.4. WISE LABs: AGN Feedback Caught in the Act? 

The correlation between the masses of supermassive black 
holes (SMBH) at the center of nearby ellipticals and their 
bulge stellar velocity dispersion remains striking in the study of 
galaxy evolution (e.g., Ferrarese & Merritt 2000). The general 
paradigm that ties these passive galaxies with their high-redshift 
progenitors is a process of merger-induced star formation which 
also fuels the SMBH (e.g., Sanders & Mirabel 1996; Hopkins 
et al. 2006). The system appears at one point as a heavily 
obscured star-forming galaxy, and later, as the SMBH accretes 
at a higher rate, a galaxy with an active nucleus. AGN- and 
starburst-induced winds eventually become strong enough to 
expel the obscuring gas and dust, briefly revealing an optical 
quasar. This short-lived “feedback” process is thought to quench 
both star formation and AGN activity, leading ultimately to a 
passive, red galaxy spheroid (e.g., Hopkins et al. 2006; Farrah 
et al. 2012). Observational evidence for this process has been 
challenging to obtain, but models provide a general picture of 
the properties of these systems. 

It is well established that the peak of this feedback activity 
occurs near z — 2 (e.g., Di Matteo et al. 2005), and that a 
multitude of galaxy types are associated with the process. SMGs 
are thought to be the IR-luminous starbursting precursors to 
systems that will evolve into a massive elliptical, and have been 
well characterized in a number of studies (e.g., Blain et al. 1999; 
Chapman et al. 2005; Magnelli et al. 2012). 

Using a three-dimensional radiative equilibrium code, 
Chakrabarti et al. (2007) predict Herschel colors for systems 
undergoing feedback from SN and AGN feedback separately 
and showed that AGN feedback alone is particularly effective 
at dispersing gas and dust. They predict that for a brief time 
(<40 Myr) while the AGN injects energy into the galaxy and 
surrounding IGM at its maximal rate, the dust is hotter than 
during the SMG starburst phase. SMGs (including the “warm” 
SMGs discovered by Herschel) typically have far-IR colors con- 
sistent with only SN feedback, while all the WLABs have Her- 
schel colors (F\\ofim/E\iOfim) at least twice as red (~0.6), 
in line with the model predictions for when AGN feedback is 
at or near its peak. While Chakrabarti et al. (2007) does not 
make predictions regarding the Lya morphology, other models 
do predict that given a clumpy gas and dust distribution, the Lya 
would also likely not be symmetrical which is consistent with 
the filamentary Lya we have observed (Cen & Zheng 2012). 
However, more detailed modeling is required. 

Finally, the surface density of WLABs is generally consis- 
tent with predicted numbers of HyLIRGs from models (e.g., 
B6thermin et al. 2011) coupled with timescale arguments for 
the feedback process and length of the SMG phase. However 
we caution that this is not a particularly strong argument given 
the uncertainty in these models as well. Nevertheless, the accu- 
mulated evidence suggesting that the WLABs are undergoing 
extremely powerful feedback is intriguing, and we will expand 
on it in future papers. 

4. CONCLUSIONS AND FUTURE WORK 

This discovery paper presents a new WISE color criterion 
with a 78% success rate in selecting z = 1. 6-4.6 dusty, mid-IR 
bright, LAEs, of which at least 37% (perhaps all) are LABs. This 
new population of rare (—0. 1 deg -2 ) galaxies, previously missed 
in narrow-area surveys, have a redshift distribution that peaks 
at z — 2.3, and total IR luminosities and dust temperatures on 
average brighter and hotter than SMGs and other high-redshift 
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dusty galaxies. This is the first systematic search technique 
to highlight bright, extremely dusty (hot) LAEs/LABs* and 
unlike optical narrow-band searches covers a large redshift 
range, the whole sky, and without contamination by low redshift 
interlopers (i.e., mistaking [O n] for Lya) making it well suited 
to providing taigets for future large spectroscopic surveys. 

The current evidence suggests that these galaxies are in a 
short-lived transition between dusty starburst and an optical 
QSO driven by the central AGN. If true, these systems offer 
a unique opportunity to investigate AGN feedback and how it 
can effect not only the galaxy but also the surrounding IGM. A 
full census of the optical spectroscopy, mid-far-IR properties, 
near-IR morphologies and catalog of the WLAEs/WLABs will 
be presented in a series of forthcoming papers. 
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